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NOTATION 


A,B,C,D,  E  ,  Abbreviations  introduced  in  Equation  (18) 
a,  b  Stretching  factors  defined  by  Equation  ( 1 5) 


c 

Ek 

f 

H 

k 

— > 
k 

L 

1,  m 
p' 

P 

Ro 

r’,  <f>,  z' 

r,  (/>,  z 

tT 

t 

u’,vy,w’ 

u,v,w 


Constant  -defined  -in  Equation  (29) 

Ekman  number 
Frequency 
Height  of  the  tank 
Summation  index 

Unit  vector  in  the  direction  of  rotation 

Radius  of  the  tank 

Eigenvalues 

Pressure 

Dimensionless  pressure 

Rossby  number 

Cylindrical  polar  coordinates 

Dimensionless  cylindrical  polar  coordinates 

Time 

Dimensionless  time 

Velocity  components  corresponding  to  r’,0,z’ 
Dimensionless  velocity  components 


ot,p  Abbreviations  introduced  in  Equation  (18) 

Abbreviation  introduced  in  Equation  (28) 

6  Ratio  H/L 


—♦ 

£ 

t 


v>  e 


Dimensionless  vorticity  vector 
Azimuthal  component  of  £ 
Stretched  coordinates  z,  r 
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60 


Kinematic  viscosity 
Density 

Abbreviations  introduced  in  Equation  (28) 
Dimensionless  stream  function 
Function  introduced  in  Equation  (28) 

Angular  velocity  of  the  tank 
Angular  velocity  of  the  cover 

=  n  -  a 

c 


Sub-  and  Superscripts: 

i,  j  Location  of  grid  point  in  the  (tj,  0) -plane 

n  Location  of  grid  point  in  time 
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ABSTRACT 


A  study  has  been  made  of  axisymmetric 
incompressible  fluid  flows  in  a  rotating  tank  when  the 
angular  speed  of  the  cover  changes  abruptly.  From 
the  initial  solid-body  rotation  a  meridional  and  an 
azimuthal  circulation  relative  to  the  moving  tank 
develop.  This  problem  is  solved  numerically  by 
means  of  a  stream  function-vorticity  formulation 
for  the  meridional  flow.  Local  fine  grids  are  used 
in  the  Ekman  and  Stewartson  layers.  No  finite  gap 
between  tank  and  cover  is  considered.  The  singular 
behavior  at  this  point  is  investigated.  The  parameters 
considered  are  the  Rossby  number,  the  Ekman 
number,  and  the  ratio  of  height  to  radius  of  the  tank. 
Temporal  and  spatial  oscillations  of  the  laminar 
flow  field  as  well  as  the  occurrence  of  cell  flows  are 
discussed. 
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1.  INTRODUCTION 

In  recent  years  the  study  of  rotating  fluids  has  found  widespread 
interest  from  both  the  theoretical  and  practical  point  of  view. 
Applications  of  this  subject  range  from  lubrication  problems  to 
centrifuge  design  to  geophysical  questions.  Physically,  rotating 
fluids  can  behave  quite  differently  froim  nonrotating  ones,  a  fact  which 
has  important  mathematical  implications.  In  the  Computation  and 
Mathematics  Department  the  program  to  simulate  viscous  fluid  flows 
by  means  of  computers  includes  the  study  of  rotating  fluids.  The 
present  analysis  is  an  outgrowth  of  this  effort. 

We  consider  the  following  laminar  flow  problem.  A  circular- 
cylindrical  tank,  completely  filled  with  an  incompressible  fluid,  is 
rotating  around  its  axis  with  a  constant  angular  velocity.  The  cover 
is  a  disk  which  can  rotate  co-axially  with  the  tank  but  at  a  different 
rate.  Initially,  both  tank  and  cover  rotate  with  the  same  angular 
velocity,  and  the  fluid  inside  behaves  as  a  solid  body.  At  a  certain 
instant,  the  angular  speed  of  the  cover  abruptly  changes  to  a  rate 
which  is  different  from  that  of  the  tank.  As  a  result,  a  meridional 
and  an  azimuthal  circulation  relative  to  the  moving  tank  develop  with 
time  and  approach  asymptotically  a  steady  state.  A  numerical 
computation  of  this  flow  is  the  subject  of  the  present  study. 

A  number  of  finite-difference  solutions,  obtained  with  the 
stream  function-vorticity  formulation  for  the  meridional  flow,  have 
recently  been  published  in  the  literature.  Steady-state  integrals  for 
flows  within  a  fixed  tank  with  a  spinning  cover  were  found  by 

1  2 

Dorfman  and  Romanenko  and  by  Pao  .  The  steady  motions  in  a 
References  are  listed  on  page  43  . 
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rotating  tank  with  a  fixed  cover  were  studied  by  Pao  and  by 
4  3 

Farris  et  al  .  Pao  also  computed  the  initial  phase  of  the  flows 
after  the  abrupt  change  of  the  cover's  rotational  speed.  Flows 

with  sources  and  sinks  on  the  boundaries  were  investigated  by 

4  5 

Farris  et  al  for  the  steady-state  case  and  by  Krause  for  the 

transient  case.  Small  deviations  from  solid-body  rotation  allow  the 

linearization  of  the  equations  of  motion.  This  approach  was  chosen 
6 

by  Rasmussen  .  The  spin-up  and  spin-down  of  a  cylindrical  tank 

'  7 

without  relative  rotation  of  the  cover  was  studied  by  Briley  and  Walls  . 

A  wealth  of  literature  exists  for  limiting  cases  which  permit 
similarity  assumptions  or  the  use  of  perturbation  methods.  The 
solutions  obtained  are  valuable  as  a  means  of  checking  the  accuracy 
of  the  numerical  techniques.  Older  papers  are  recorded  in 

8  9 

Schlichting's  book  ,  newer  papers  in  Greenspan’s  monograph  .  In 
particular,  attention  is  called  to  the  paper  by  Pearson***.  The  spin-up 

9 

and  spin-down  problems  were  studied  by  Greenspan  and  by 
Euteneuer  et  al**.  Details  of  some  of  the  papers  are  discussed  in 
context  with  our  results. 

The  problem  outlined  in  the  beginning  of  this  section  is  solved 
by  means  of  a  stream  function-vorticity  formulation.  The  numerical 
technique  is  essentially  that  presented  by  the  authors  in  Reference  12. 
Grid  systems  of  41x41  =  1681  and  51x51  =  2601  mesh  points  with 
unequal  spacing  represent  the  flow  field.  Thus,  local  fine  grids  can 
be  used  in  regions  with  high  vorticity  gradients.  The  parameters  of 
the  problem  are  the  Rossby  number  Ro,  Ekman  number  Ek,  and 
the  aspect  ratio  6  of  the  tank. 

The  assignment  of  values  to  these  parameters  is  guided  by  the 
objective  of  studying  flows  in  a  tank  with  nonvanishing  angular  speed. 

It  is  of  advantage  to  solve  the  equations  of  motion  in  a  rotating  frame 
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when  studying  the  elliptic  or  hyperbolic  flow  behavior  in  space  or  in 
time  for  small  Rossby  numbers.  This  can  be  seen  immediately  for 
the  linearized  case,  Ro  =  0.  Then  the  nondimensional  vorticity 
transport  equation  takes  the  form 

(  — y-  -  2Ek  v2  +  Ek2  v4)  v2  C+(2k*  v)2  £=  0  , 

St2  9t 

— + 

where  k  is  the  unit  vector  in  the  direction  of  rotation.  The  other 

quantities  are  described  in  Section  2.  When  Ek  is  small  the  hyperbolic 

form  of  the  equation  in  time  is  revealed.  For  steady  motions  we 

g 

arrive  at  the  Taylor-Proudman  theorem  if  Ek  =  0  , 


2.  FORMULATION  OF  THE  INITIAL -BOUNDARY 
VALUE  PROBLEM 


We  assume  a  laminar  axisymmetric  flow  of  an  incompressible 
fluid  in  a  circular- cylindrical  tank  of  radius  L  and  height  H,  which 
is  spinning  with  constant  angular  velocity  Q .  At  time  t'=0  the 
cover  impulsively  starts  to  rotate  with  a  different  but  constant 
angular  velocity  Oc.  Cylindrical  polar  coordinates  r',  </>,  z’  are 
used  with  the  corresponding  velocity  components  u’,  v’,  w*  in  a 
reference  frame  rotating  with  the  tank.  Under  the  restriction  of 
axisymmetry,  a/&0  =0,  the  Navier-Stokes  equations  and  the  equation 
of  continuity  are 


uV+U,U,r,  +  W,U,z’_V?(20  +  7r^‘  ^p,r,+  ^UVr’+(7V,+U,zfz,]’  (1) 

vV+u’V,r,  +  w’v,z,  +  u,(2Q+rr)  =  t/[v,rv+(P')  .  +  v’zV]’  (2) 
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wV.  +  u’w'  ,  +  w'w' 
t’  r*  .  s 


i  w  , 

.pPV^^rV^^’zvl' 


(3) 


(4) 


Here,  p',  p,  and  v  are  the  pressure,  the  constant  density,  and  the 
constant  kinematic  viscosity,  respectively.  Prior  to  the  sudden 
change  of  rotation  of  the  cover,  the  entire  fluid  is  at  rest  relative  to 
the  spinning  tank:  . 

tf  <0:  u’ s  0,  v’  h  0,  w'  *  0  .  (5) 

After  the  change  of  rotation,  t'*  0,  the  boundary  conditions  are 

z'  =  0,  0  sr'  sL:  u’  =  0,  v'  =  0,  w'  =  0, 

z'  =  H,0sr'sL:  u'  -  0,  v'  =  -cor',  w'  =  0,  (6) 

r’  =  L,  0  < z*  < H:  u'  =  0,  v'  =  0,  w'  =  0  , 

where  co=  O  -  O  .  It  is  convenient  to  introduce  the  following 
dimensionless  variables; 

t'=t /O,  r'=Lr,  z'=Hz,  (u',  v’)=coL(u,v),  w'=coHw,  p'  =  pcoOL  p  (7) 
and  the  characteristic  numbers 


Ro  = 

CO 

O 

(Ross by  number), 

Ek  = 

V 

(Ekman  number), 

(8) 

OH 1 

'  6  = 

H 

L 

(aspect  ratio), 

The  axisymmetry  of  the  motion  permits  the  stream  function-vorticity 
formulation  of  the  meridional  flow,  If  $  designates  the  dimensionless 
stream  function  and  C  the  azimuthal  component  of  the  dimensionless 
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vorticity  veqtor  £  ,  where 

u  =  i  0Z  ,  w  =  -  i  «/>r  ,  (9) 

£=uz-62wr,  (10) 

Equations  (1)  through  (4)  are  reduced  to 

^+Ro[<ut)  +(wj) 0  + tzz]  .  (11) 

jr  r 

v  +Ro[(uv)  +(wv)  +2^]+2u  -  Ek[62(vrr  +  ivr-^)+v  ]  ,  (12) 

r 

|[«2(*rr-i*r)  +  *zz]  =  t.  (13) 

The  corresponding  boundary  conditions  for  t  *  0  are 

z  =  0,  0sr<  1:  0  =  0,  0  =  0,  v  =  0, 

z 

z  e  1,  0  s  r  s  1:  0=0,  0  =  0,  v  =  ~r, 

Z  (14) 

r  =  1,  0  £  z  <  1:  0=  0,  0r  -  0,  v  -  0, 

r  =  0,  0  s  z  <  1:  0=0,  £  =  0,  v  =  0. 

The  last  conditions  for  the  values  on  the  centerline  follow  from  the 
axisymmetry  of  the  flow,  which  allows  restriction  to  half  the 
meridional  plane. 

At  the  corner  r  =  1,  z  =  1  a  discontinuity  in  v  occurs  which 

causes  the  shear  stresses  lim  yp(|^)  and  lim  yp(|^-)  to 

r-4  “z  z=l  z-*l  °r  r=l 

be  unbounded  at  that  point.  In  reality,  a  small  gap  between  cover  and 

container  always  exists.  Thus,  this  singularity  is  avoided.  The 

13 

matter  was  discussed  by  Schmieden  for  slow  motion  with  no 
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meridional  flow.  If  the  gap  is  taken  to  be  infinitesimal  (as  in  this 
work),  the  torques  exerted  on  the  cover  and  the  tank  are  always 
logarithmically  singular.  Although  the  singularity  becomes  more  and 
more  localized  with  decreasing  Ekman  number,  the  torques  are  always 
infinite.  The  implications  for  the  finite-difference  scheme  are 
discussed  below. 


3.  NUMERICAL  INTEGRATION 


For  finite-difference  techniques  it  is  of  advantage  to  lay  out  a 
network  which  is  dense  in  regions  of  high  vorticity  gradients.  Such 
regions  usually  occur  as  boundary  layers  near  solid  walls.  For  the 
present  problem  the  following  coordinate  transformation  is  introduced: 


r  =  0  +  a  sin  7T  0  , 
z  =  n  -  b  sin  2t r  r]  , 


(15) 


where  the  stretching  factors  a  and  b  are  chosen  from  the  intervals 
0  =sa  s0.26  and  0. 10  £  b  s;  0. 13.  The  region  of  integration  is 
represented  by  a  grid  with  mesh  points  at  17=  =  (i- 1) A 77 , 

0=  0-  =  (j-1)  A  0  ,  where  i,  j  =  1, . . , ,  41  or  51.  (See  Figure  1.) 

J 


After  the  equations  of  motion  are  transformed  according  to 
Equation  (15),  they  are  solved  with  the  following  finite-difference 
scheme.  The  linear  differential  operators  of  Equations  (11)  and 
(12),  except  for  the  Coriolis  terms,  are  replaced  by  the  Dufort- 
Frankel  approximation.  The  term  -2vz  and  the  nonlinear  operators 
are  expressed  by  central-difference  formulae.  Then,  Equations  (11) 
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Figure  1  -  Grid  System  with  41  x  41  Mesh  Points 
and  Stretching  Factors  a  *  0.2,  b  =  0.1 
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and  (12)  yield,  for  the  (n+l)^  time  step, 

(1  -iE^tC^  ^"^C^  +  aAt  [Ro(a .  3 1'1-1*  J-^i,  3+lCi,  3+1 

n  „n  n  '  n  n  n  n 

w.  i  .C.  ,  • -w. -  .  v.  .  v. . ,  .  -v.  *  . 

+  8  Jlhrthl _ 1+13^+13.  +  i>  3  8  1+1 3  thl) 

pi  2An  t»  Pi  An  ' 

,  y 

vn  . -vn  . 

x  A  i+1  3 _ H_3  .  rn  /A  >»n  x  n  )»n  X  n  -0-1 


LJllL+  E,  (A.C«  ,^1  +  B.c  ii+  C,  .c 
An  kv  3  l,  j+l  3  iJ+1  i,3  i,i 


VM1=(1  -2EkAtCi,j)'1  t\l2^[Ro(«.  Ui-+1-TVA"e1  Ui’i+lVl’+1 


n  n  n  n 
w.  -  .v.  i  .-w.  *  .v.  . 
i-l.j  i-l,3  i+1,3  1+ 


2An 


2  n  n  x 
—  u.  .v.  .)  -  2u.  .. 
i>3  i,3  1,3 


+ Wi,  J+i+Bjv"  i-i+ci,  ivij1+Divi+i,i+Eivi-i,  j)  n  •  (i7) 


Here,  the  following  abbreviations  are  used: 


a,  = 


,  ^4  S 


j  1+  air  cos  7T  0 .  ’  Ki  ~  1- 2ir  bcos  2ir  n^  * 


r.  s  e.  +  asin7r  6.  , 

3  3  3  ’ 

s  rt  2 ,2  3  ,  .2 

Oa .  2  it  o  a,  a  sinir  0,  o  a . 

A  s  (  — -3-)  +  - J - J-  +  — - 3 — 

j  K  &6  '  2A0  2r  A0 

J 


»  2.2  3  .  ,2 

Oa .  2  7T  o  a .  a  sinff  0.  o  «. 

•p  ss  / _ 3  \  _ _ 3  _ 1.  _ _ 3_ 

3  V  A0  1  2A0  2r.A  0  ’ 

J 
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6a.  2 

-  (  — 1 ) 

1  A0  ' 


&  2 

’ 


j8.  2  2772  b  sin  2  77  r). 

(_J_)  .  - I - L 

'  Atj  ■  At? 

0.  2  2t72/3  2  b  sin  2t7  17. 

(_JL  )  +  - i - L 

v  A  n  at? 


The  velocity  components  of  the  meridional  flow  are  obtained  from 


u  -  ^i+l,3~^i-lj 

ij  3  r.  2  A?? 


-  “i  ^i,j+l~^ij-l 


■or  = - —  7  - - 

’  fcj  r.  2a  0 


The  Poisson-type  Equation  (13)  is  approximated  by  the  five  - 

point  formula  which  yields,  for  j/>.  .  , 

3 

1  6  a  •  2  j8.  2  j 

lb.  .=i  [(-—-)  +  (tM  1’  rD.  ib.  ,  .  +  E.  it .  *  . 

2  L  Ar?  At?  J  L  iu+l,j  1  ^ i-l,  J 

6| .  62a . 

+  ^Aj~  r^Ae^i,  j+l+^Bj+  r.A0  ^i,j-l  j  ^  ^ 

This  system  of  algebraic  equations  is  solved  with  Gauss -Seidel  line 

overrelaxation  applied  along  lines  of  constant  0  .  The  overrelaxation 

til 

factor  is  1.  78.  The  iteration  is  halted  after  the  k  iteration  if,  at 


each  grid  point. 


-  €) 1 K 10’4  'C/ 


At  the  solid  boundaries  a  one-sided  first  order  difference  equation  is 
used  to  compute  the  vorticity  £ .  Two  such  equations  were  tried. 
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the  simple  one  which  is  used  by  most  authors  (at  the  bottom  of  the 
tank,  for  instance), 


A 


’1,3 


1 12 


2,3 


(22) 


12 

and  another,  which  was  found  in  a  study  involving  curved  boundaries 
to  be  superior  with  regard  to  numerical  stability, 


cl,j  =  4rj  (  tr,  )  (*2,j  +  4*’3,j  ■  *4,J) 


(23) 


In  this  paper  Equation  (23)  is  used  although  it  showed  no  advantage 
over  Equation  (22). 

It  can  be  seen  from  Equations  (16)  through  (20)  and  Equation 
(23)  that  the  computations  at  the  inner  points  do  not  require  knowledge 
of  quantities  at  the  singular  point.  (The  subscript  i+1,  j+1  does  not 
appear  in  these  equations.) 

The  integration  process  is  carried  out  in  the  following  way:  The 

vorticity  and  subsequently  the  azimuthal  velocity  v!14^  are 

computed  at  the  inner  points  according  to  Equations  (16)  and  (17).  . 

n+1 

The  calculation  of  if).  .  follows  with  the  aid  of  Equation  (20).  The 

l,  3 


cycle  then  concludes  with  the  calculation  of  the  surface  vorticity. 

The  maximum  stable  time  step,  Atmax  ,  beyond  which  numerical 
instability  occurs,  is  determined  by  increasing  the  time  step  until 
oscillations  from  point  to  point  in  the  £  -values  appear. 

The  accuracy  of  the  computations  is  checked  by  using  different 
space  increments.  This  is  easily  done  by  varying  the  stretching 
parameters  a  and  b  defined  in  Equation  (15).  It  is  found  for  the 
two  cases  Ro  =  1,  6=  1,  Ek  =  0.01  and  0.001  that  solutions  obtained 
with  different  stretching  values  agree  well  except  near  the  singular 
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point.  The  influence  of  this  singularity  on  the  solution  is  discussed 
in  Section  4.  Another  way  to  change  the  space  increments  is  to 
vary  the  number  of  grid  points.  For  Ro  =1,  Ek  =  0.  001  two  grid 
systems,  41  x  41  and  51  x  51,  are  used.  The  agreement  is  good. 


4.  RESULTS 


As  already  noted,  our  main  interest  is  focused  on  fluid  motions 

in  a  rotating  tank,  that  is  flows  with  Ro  <  »  .  The  examples  selected 

for  computation  are  compiled  in  Table  1.  All  cases  are  started 

from  solid-body  rotation  at  t  =  0  and  are  continued  to  an  almost  steady 

state  at  t_.  ,  ,  where  t_,.  ,  is  the  earliest  time  at  which 

Final  ’  Final 


|j/)n  .  -  ^)n_1  (<  10~4  l^”1!  (24) 

is  satisfied  throughout  the  field.  The  computations  were  performed 
in  double  precision  on  an  IBM  360-91  computer  and  in  single 
precision  on  a  CDC  6700  computer.  Pictures  of  the  flow  field  were 
made  with  a  Stromberg-  Carlson  SC -4020  charactron  plotter.  / 

a.  The  Almost  Steady  Case 

Two  cases  may  be  distinguished  if  the  nonlinear  inertial  terms 
are  neglibible.  For  large  Ekman  numbers  the  pressure  force  is 
essentially  balanced  only  by  the  friction  force  (slow  motion).  For 
Ro  «  Ek  «  1  a  balance  is  maintained  among  the  Coriolis,  friction, 
and  pressure  forces. 

14 

Slow-motion  solutions  have  been  obtained  by  Hort  for  Ro  =  1, 

3 

6=  1,  Ek  =  *  .  Pao  checked  numerically  for  Ro  =  1,  6=1  that  the 
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TABLE  1 

Compilation  of  the  Calculated  Examples 


No. 

Ro 

Ek 

6 

a 

b 

GRID 

1 

10-5 

0.01 

1 

0.2 

0.1 

41  x  41 

2 

10“5 

0.001 

1 

0.2 

0.1 

3 

1 

O 

H 

0.01 

1 

0.2 

0.1 

4 

1 

0.01 

1 

0.0 

0.1 

5 

1 

0.01 

1 

0.26 

0.13 

6 

1 

0.005 

1 

0.15 

0.1 

7 

1 

0.0025 

1 

0.2 

0.1 

8 

1 

0.00125 

1 

0.2 

0.1 

9 

1 

0.001 

1 

0.2 

0,1 

10 

1 

0.001 

1 

0,26 

0,13 

11 

1 

0,001 

1 

0.0 

0.0 

51  x  51 

12 

1 

0.0002 

1 

0.2 

0.1 

51  x  51 

13 

1 

0.01 

3 

0,0 

0.1 

41  x  41 

14 

1 

0.001 

3 

0.2 

0.1 

15 

1 

0.01 

1/2 

0.2 

0.1 

16 

1 

0.01 

1/3 

0,2 

0.1 

17 

4 

0.01 

1 

0.2 

0.1 

18 

4 

0.004 

1 

0.2 

0.1 

19 

10 

0.01 

1 

0.2 

0.1 

13 


slow-motion  approximation  is  applicable  for  0. 125  <  Ek  £  »  . 

The  linear  theory  predicts  for  the  case  Ro  «  Ek  «  1  three 

distinct  regions:  The  Ekman  layers  at  the  cover  and  at  the  bottom  of 

the  tank,  the  Stewartson  layer  at  the  side  wall,  and  the  geostrophic 

interior,  for  which  the  Taylor-Proudman  theorem  of  inviscid  fluids 
16 

holds  .  Numerical  results  verify  this  notion  and  reveal  its 

limitation.  Figure  2  shows  lines  of  constant  ip ,  £  ,  andv-values  for 

Ro  =  10“4,  Ek  =  0.01,  6=1  and  for  Ro  =  10“5,  Ek  =  0.001,  6  =  1. 

The  increments  of  the  ip,  C  >  and  v -values  are  recorded  in  Table  2. 

The  case  Ro  =  10“^,  Ek  =  0.001,  in  particular  the  lines  of  constant  £, 

illustrates  clearly  the  existence  of  the  three  regions  predicted  by  the 

-5 

linear  theory.  Additional  results  obtained  for  Ro  =  10  ,  Ek  =  0.01, 

6=1  agree  so  well  with  the  case  Ro  =  10"^,  Ek  =  0.01,  6=1  that 

the  patterns  are  indistinguishable.  The  increase  of  the  Ekman 

number  from  0.  001  to  0.01  shows  that  the  (almost)  inviscid  interior 

ha6  vanished.  If  we  neglect  the  influence  of  the  side  wall,  an  analytic 

solution  from  the  linear  theory  is  available.  This  is  the  well-known 

Ekman  solution.  In  Figure  3  the  numerical  values  for  v  at  the 

centerline  r  =  0  are  compared  with  those  of  the  Ekman  solution. 

Again,  the  analytic  values  agree  well  with  the  numerical  output  for 

Ek  =  0. 001,  whereas  for  Ek  =  0. 01  the  numerical  data  reveal  the 

influence  of  the  side  wall.  A  distinct  asymmetry  is  displayed  between 

cover  and  bottom  for  Ek  =  0. 01. 

The  Stewartson  layer  at  the  side  wall  is  studied  with  a 

g 

perturbation  method  and  compared  with  the  numerical  output.  Two 

different  layers,  one  inside  the  other,  must  be  distinguished:  a 

1/4  1/3 

layer  of  thickness  Ek  ,  and  a  layer  of  thickness  Ek  inside  the 

first  layer  and  adjacent  to  the  wall.  The  Ek*/4  -layer  is  represented 
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TABLE  2 

Data  for  the  0 ,  C  >  and  v  Patterns  in  the  Figures 


No. 

Figure 

Sign  of  0 

Near  r=l,  z-1 

A0 

Sign  of  C 
on  the  Cover 

AC 

Av 

2 

2 

- 

0.001 

- 

1.00 

0.05 

3 

2 

- 

0.001 

- 

0.25 

0.05 

5 

6 

- 

0.001 

- 

0.25 

0.05 

7 

6 

- 

0.001 

- 

0.25 

0.05 

8 

7 

- 

0.001 

- 

0.50 

0.05 

10 

7 

- 

0.001 

- 

0.50 

0.05 

12 

8 

- 

0.001 

- 

1.00 

0.05 

19 

8 

+ 

0.001 

+ 

0.50 

0.05 

17 

9 

+ 

0.001 

+ 

0.25 

0.05 

18 

9 

+ 

0.001 

+ 

0.25 

0.05 

13 

10 

- 

0.0001 

- 

0.25 

0.05 

16 

10 

- 

0.001 

- 

0.25 

0.05 

14 

11 

- 

0.001 

- 

0.50 

0.05 

2 

15 

- 

0.001 

- 

-- 

0.05 

10 

19 

- 

0.001 

- 

-- 

0.05 

12 

20 

- 

0.001 

- 

-- 

12 

21 

- 

0.001 

- 

-- 

0.05 

19 

22 

+ 

0.001 

+ 

0.50 

0.05 

The  stream  function  0  is  specified  to  be  zero  at  the  boundaries.  The 
vorticity  C  is  zero  on  the  centerline.  The  azimuthal  velocity  v 
is  zero  at  all  boundaries  except  at  the  cover,  where  v  =  -r. 
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Figure  3  -  Comparison  of  the  Ekman  Solution  with 
the  Numerical  Results.  -(v/r)  q  is 
Plotted  versus  z. 
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by  a  closed-form  solution,  whereas  the  Ek1' J -layer  must  be 
expressed  by  a  series  expansion.  For  simplicity,  the  velocity 
components  are  presented  for  6=1: 


u  =  |  Ek1/2  e"^2  °  -  \  Ek1/2  L  (xk)TT  cos  k7r  z  , 

k- 1 


r  X  - /2  a  1/6  $  , 

v  =  -  o“  o’  e  +  Ek  '  Z  Xl.  cos  kir  z  , 

22  k=l  k 


w  =  _  1  Ek1/2  .  i  Ek1/4  (z  -  i )  e‘/2  ° 


1  /6  00 

+  Ek  '  £  X.  sin  kff  z  , 

k=l  K 


where 


z  1xk  2  .  ,/3  ,  "yk^2  T 

jfkT  yksin(TykT>e 

a  =  (1  -  r)/EkX/4  ,  T  =  (1  -  rJ/Ek1/3  , 


7k  =  (2k  7r)' 


The  first  term  in  v  and  w  is  the  geostrophic  mode  of  the  interior. 

In  Figures  4  and  5  the  v-  and  w-components  are  plotted  against  r. 

The  numerical  results  are  compared  with  the  analytic  data  for 

1/4 

a)  the  geostrophic  mode  and  the  Ek  '  -law,  and  b)  the  geostrophic 
1/4 

mode,  the  Ek  '  -law,  and  the  first  two  nonzero  terms  of  the 
1/3 

Ek  '  -series.  For  the  v-component  both  analytic  curves  are  in  good 

agreement  with  the  numerical  data.  For  the  w-component  the 
1/4 

Ek  -law  alone  is  insufficient  to  describe  the  Stewartson  layer.  In 
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Ml] 


NUMERICAL,  Ro=10  s, 

Ek-10~3 

ANALYTIC,  6E0STR0PHIC 
INTERIOR  +  Ek'/' 

ANALYTIC,  6E0STR0PHIC 
INTERIOR  +  Ek'/4  +  Ek* 
(FIRST  TWO  NON— VANISHIN6 

TERMS) 


Figure  4  -  Analytic  and  Numerical  Results  for  the  Stewartsoii 
Layer,  -v  is  Plotted  Versus  r  at  z  =  1/2 
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NUMERICAL,  »»-1l5,  Ek-1»' 
ANALYTIC,  6E0STI0RHIC  INTERIOR  +  Ek 
ANALYTIC,  6E0STR0RNIC  INTERIOR  +  Ek 
+  Ek'7’  (FIRST  TWO  1011- 
VANISHIN6  TEINS) 


Figure  5  -  Analytic  and  Numerical  Results  for  the 
Stewartson  Layer,  w  is  Plotted 
Versus  r  at  z  *  1/2 

20 


1  1/4 

fact,  at  z  =  g  »  the  contribution  of  the  Ek  -law  vanishes. 

The  inclusion  of  the  nonlinear  terms  into  the  equations  of 

motion  modifies  the  flow  characteristics.  Figures  6  through  8 

show  flow  patterns  for  Ro  =  1,  6  =  1,  and  various  Ek.  The  output 

3 

for  Ek  =  0.01  can  be  compared  with  Pao’s  computation  .  The 

overall  streamline  patterns  agree  well  with  each  other,  although 

deviations  occur  near  the  singular  point.  This  problem  is  discussed 

in  more  detail  below.  In  Figures  8  and  9  flow  patterns  are  displayed 

for  Ro  =  10  and  4,  6  =  1,  and  two  different  Ekman  numbers.  Cell- 

type  motions  occur  which  are  expected  around  Ro  =  4  according  to 

17 

similarity  solutions  .  Results  in  Figures  10  through  12  exhibit  the 

effect  of  various  6  ,  In  the  limit  6  -*  0  the  solutions  approach  the 

15 

similarity  solutions  obtained  by  Lance  and  Rogers  ,  This  is 
demonstrated  in  Figure  12  by  a  graph  of  -v/r  versus  z  at  r  =  0. 

The  results  for  different  Ro,  Ek,  6  have  one  common  feature. 
Below  a  certain  value  of  the  Ekman  number,  spatial  undulations  of 
the  streamlines  occur  indicating  the  transition  from  an  elliptic -type 
to  a  hyperbolic -type  solution.  This  distinction  is  most  apparent  in 
the  two  linear  cases  discussed  at  the  beginning  of  this  section.  The 
slow- motion  solution  is  clearly  elliptic,  whereas  the  Ekman  solution 
and  the  geostrophic  mode  are  parabolic  (as  an  asymptotic  limit  of 
the  hyperbolic  time -dependent  flow).  For  Ro  =  1,  6=1  the  critical 
Ekman  number,  at  which  the  transition  occurs,  is  0.005.  This 
value  decreases  for  Ro  >  1  and  increases  for  Ro  <  1.  Variation  of 
6  appears  to  have  only  a  minor  effect  on  the  value  of  the  critical 
Ekman  number. 

In  figure  13  a  pathline  is  displayed  for  Ro  =  1,  Ek  =  0. 0002, 

6  =  1  in  thfe  rotating  frame.  The  starting  point  is  arbitrarily  chosen 
at  r  =  0. 4&9,  z  =  0. 959.  The  fluid  particle  follows  a  trajectory  which 
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22 


23 


24 


Figure  10  -  Lines  of  Constant  i}i,  5,  and  v  for  Ro  = 
Ek  =  0.01,  6  =  3  and  6  =  1/3  at  almost  Steady  State 


V 


Figure  11  -  Lines  of  Constant  i|/,  and  v  for  Ro  =  1,  Ek  =  0.001,  6=3 

at  almost  Steady  State 
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0 


1.0 


z 


Figure  12  -  -(v/r)  q  as  a  Function  of  z  for  Various  6 
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Figure  13  -  Pathline  for  Ro  =  1,  Ek  =  0.0002,  5=1.  Perspective  View  at  almost 

Steady  State 


consists  of  a  downward  spiral  with  almost  no  radial  variation  and 
an  upward  spiral  near  the  side  wall  with  only  one  revolution.  The 
time  for  the  particle  to  return  to  the  vicinity  of  its  initial  position 
is  about  22  revolutions  of  the  tank  (t  =  140). 

The  singularity  of  the  flow  at  r  =  1,  z  =  1  requires  special 

attention.  In  Figure  14  the  function  -  —  (•—)  is  plotted  against 

r  az  Z=1 

(1-r)  on  a  double -logarithmic  scale  for  various  values  of  Ro,  Ek,  and 
6  .  The  dimensional  quantity  a  v’/a  z'  becomes  unbounded  at  the 
singular  point  (r ’  =  L,  z*  -  H)  according  to 


lim 

r'-L 


-  (|^)  ,  =  c 

az  z  =H 


Leo 

L-r’  ’ 


c  0.  73  , 


(29) 


where  c  does  not  depend  on  Ro,  Ek,  or  6 .  Hence,  the  singularity 
is  a  local  effect  which  can  be  considered  in  a  numerical  scheme  by 
building  in  a  local  series  expansion.  On  the  other  hand,  the  torques 
exerted  on  the  cover  and  the  wall  are  always  logarithmically  singular. 
Figure  14  reveals  that  for  Ro  =  1,  Ek  =  0.01  the  grid  system  with 
a  =  0,  b  =  0. 1  gives  a  poor  representation  of  the  near-corner  region. 

b.  The  Transient  Case 

Unsteady  rotating  flows  exhibit  parabolic  or  hyperbolic 

14 

properties  in  time.  The  slow-motion  solution  ,  for  example,  is 
a  pure  diffusion  process  and  is,  thus,  parabolic.  The  linear  theory 
of  inviscid  flow,  in  which  the  pressure  gradient,  the  Coriolis  force, 
and  the  local  acceleration  balance  each  other,  reveals  flows  of 

9 

hyperbolic  nature  .  This  can  be  seen  immediately  for  the  vorticity 
equation  cited  in  the  introduction.  If  Ek  =  0,  an  infinite  but  countable 
number  of  modes  can  be  obtained.  Their  frequencies  are 
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0.01  0.02  0.04  0.1  0.2  0.4  1 


0-0 

0  Ro  =  10“5  Ek  =  10^  6=  1 
@  EKMAN  SOLUTION 
®  Ro  =  1.  Ek  =  0.01  6=1.  a  =  0,  b  =  0.1 

0  Ro  =  1,  Ek  =  0.01,  6=1,  a  =0.26.  b  =  0.13 

©  Ro  =  1,  Ek  =  10"3,  6=1 

©  Ro  =  1.  Ek  =  0.01,  6=  1/3 

©  Ro  =  1,  Ek  =  0.01,  6=3 

Figure  14  -  -  -^-(3v/3z)  as  a  Function  of  1-r  for  Various 

Z  Ro,  Ek,  and  6 
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(30) 


f  -  2m 

J~2  2 

Y  -t  +  m 

where  l  and  m  are  the  eigenvalues  for  the  r-  and  z-components, 

respectively.  The  waves  connected  with  these  modes  are  called 
g 

’’inertial  waves”  .  The  zero  mode,  which  is  time-independent, 
represents  the  geos  trophic  motion. 

Modifications  of  the  inviscid  flows  occur  in  the  Ekman  and 
Stewartson  layers  in  order  to  satisfy  the  boundary  condition  of 

-5 

no-slip.  In  Figure  15  a  sequence  of  ip  -  and  v-patterns  for  Ro  =  10  , 

Ek  -  0.001,  6=1  over  the  transient  period  is  presented.  A  more 
detailed  time  history  for  the  same  case  is  recorded  in  Figure  16, 
where  -(C/r)  q  at  the  cover  and  at  the  bottom  is  plotted  against 
time.  After  the  sudden  change  of  the  cover's  angular  speed,  a 
boundary  layer  on  the  cover  develops  whose  -(C/r)  q  -value  reaches 
a  peak  at  t  =  1.5.  This  time  agrees  with  the  spin -up  theory  which 

9 

estimates  t  1  for  the  initial  phase  .  Afterwards,  inertial 
oscillations  are  visible.  They  appear  in  the  streamline  patterns  in 
the  form  of  temporal  oscillations  and  produce  cell-type  motions  in 
the  center  of  the  tank  (Figure  15,  t  =  3.  98,  9.  98).  With  increasing 
time  the  inertial  oscillations  are  damped,  the  local  cell  vanishes, 
anf  the  flow  reaches  a  steady  state.  Both  curves  in  Figure  16 
approach  asymptotically  the  value  A  Ek"*^  =  15.  811  of  the  Ekman- 

a 

layer  solution. 

With  increasing  Ekman  number  the  influence  of  viscosity  is 
felt  in  the  interior.  In  Figure  17  the  function  -(C/r)  q  at  the 
cover  and  at  the  bottom  is  plotted  against  time  for  Ro  =  10“®  , 

Ek  =  0.01,  6=1.  Inertial  oscillations  have  almost  vanished.  In 
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Figure  16  -(?/r)r=Q  at  Bottom  and  Cover  as  a  Function  of  t  for  Ro  =  lcf5,  Ek  =  0.001 


Figure  18  the  flow  behavior  away  from  the  boundary  layers  is 
displayedo  The  function  -(v/r)r_Q  at  z  =  1  is  plotted  against  time  for 


both  Ek  “■  0.001  and  Ek  "  QoOl  (Ro  =  10  6“  1).  Inertial  oscillations 

are  observed  only  for  the  flow  with  the  smaller  Ekman  number„  No 

single  mode  can  be  identified.  This  indicates  the  existence  of  nonzero 

eigenvalues  l  in  Equation  (30). 

As  long  as  the  Ekman  number  is  sufficiently  small,  that  is, 

Ek  «  0.01,  inertial  oscillations  also  occur  if  nonlinear  effects  are 

present.,  For  Ro  ™  1,  Ek  =  Q.Q01,  6  =  1  a  time  sequence  of  \p-  and 

v -patterns  is  shown  in  Figure  19.  Again,  a  cell  is  visible  in  the  center 

of  the  tank  as  in  the  linear  case  Ro  =  10"°  ,  Ek  -  0.001.  A  new 

phenomenon  in  Figure  19  is  the  appearance  of  a  local  region  of 
3 

positive  V.  Pao  computed  the  case  Ro  -  1,  Ek  --  0.  001  up  to  t  *  4. 2. 
His  streamlines  agree  well  with  ours,  but  his  picture  does  not  show 
the  local  cell  at  t~4. 2. 

For  Ro  ~  1,  Ek  -  0.001,  6=1  a  computer-generated  movie 
has  been  made  which  shows  the  transient  stage  for  ip  and  v.  The  movie 
clearly  reveals  the  time  oscillations  which  cannot  easily  be  detected 
in  Figure  19. 

Higher  order  modes  are  observed  for  decreasing  Ekman 

number.  In  Figures  20  and  21  the  0-  and  v -patterns  are  displayed  for 

the  transient  period  of  the  case  Ro  -  1,  Ek  =  Q.0002;ij  6  =  1.  At  t  =  4.  50 

and  11.25  two  cells  at  the  centerline  are  visible.  The  whole  sequence 

shows  three  distinct  time  periods  when  local  cells  are  present.  It 

3 

may  be  mentioned  that  Pao  computed  this  case  up  to  t  =  3.0.  His 
strong  undulations  of  the  streamlines  near  the  singular  point  are  not 
verified  by  our  calculations. 
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as  a  Function  of  t  for  Ro  =  10  ,  <5  =  1,  Ek  =  0.001  and  0.01 


Figure  19  -  Lines  of  Constant  ^  and  v  for  Ro  -  1,  Ek  -  0.001,  5  =  1  at  Various  t.  The  almost 
Steady  State  at  tFINAL  =  102  is  Included  in  Figure  7 


Figure  21  -  Lines  of  Constant  v  for  Ro  =  1,  Ek  =0.0002,  6  =  1  at  Various  t.  The  almost 

Steady  State  at  =  245  is  Included  in  Figure  8 


Finally,  a  time  sequence  of  constant  ip ,  £  ,  and  v  curves 
is  presented  for  Ro  =  10,  Ek  =  0.  01,  6  -  1  in  Figure  22.  At  this 
high  Rossby  number  the  oscillations  are  shifted  toward  the  sidewall. 


5.  CONCLUSIONS 

Simple  geometrical  tank  configurations  with  simple  initial  and 
boundary  conditions  for  the  fluid  inside  can  generate  complicated  flow 
circulations.  In  the  transient  period  inertial  waves  are  visible  in  the 
form  of  oscillating  streamlines  and  cell  motions  without  preferred 
mode.  In  the  steady-state  case  monotonic  and  undulating  streamline 
patterns  are  distinguished. 

The  numerical  analysis  which  was  developed  for  nonrotating 
motion  also  works  well  for  rotating  flows  with  hyperbolic  features.  The 
numerical  calculations  have  been  restricted  to  flows  with  Ek  ^  0.0002 
(for  Ro  =  1)  since  the  assumption  of  axisymmetry  does  not  seem  to  be 
justified  for  smaller  Ek. 

The  linear  theory  developed  in  literature  on  the  basis  of 
perturbation  methods  is  restricted  for  the  case  Ro  «  Ek  «  1  to 

_3 

Ek  js  10  (for  6=1)  and  for  slow  motion  according  to  Reference  3  to 
Ek  >  0. 125  (for  Ro  =  1,  6=1). 

As  6  -»  0  the  computed  values  approach  the  similarity  solutions 
for  two  infinite  disks. 

The  infinitely  thin  gap  between  side  wall  and  cover  causes  a  flow 
singularity  which,  although  local  and  very  weak,  results  always  in 
an  infinite  torque.  For  practical  applications  a  nonvanishing  gap  width 
must  be  considered.  The  strength  of  the  torque  appears  to  depend 
crucially  on  this  width  (see  also  Reference  13). 


41 


42 


Steady  State  at  t  -  50  is  Included  in  Figure 
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